Bond and Neel order and fractionalization in ground states of easy-plane 

antiferromagnets in two dimensions 
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We describe the zero temperature phases and phase boundaries of an effective lattice model for 
frustrated, easy-plane, spin S — 1/2 quantum antiferromagnets. The model contains states with Neel 
and bond-centered charge order and fractionalization in different limiting regimes, and we describe 
how these orders compete at intermediate coupling. The complex Berry phase terms render this 
model unsuitable for direct simulation; however, a duality mapping leads to a model of interacting 
current loops with positive weights, which is then studied by Monte Carlo simulations. 



Understanding the zero temperature phase diagram of 
frustrated, spin S = 1/2, quantum antiferromagnets in 
two dimensions Q] is one of the central problems in the 
theory of correlated electron systems. Such an under- 
standing is expected to prove valuable in a description 
of the cuprate superconductors. Moreover, a number of 
frustrated two-dimensional antiferromagnets with a non- 
magnetic ground state have been discovered^, ||]: these 
could possibly be doped with charge carriers in the fu- 
ture, leading to exciting possibilities for new physics. 

Some understanding of the different classes of possible 
ground states of insulating antiferromagnets has emerged 
in recent years. States can be classified by their posses- 
sion of one or more distinct types of 'order' Q: mag- 
netic order involving breaking of spin-rotation invariance 
(as in a Neel state), bond order arising from a sponta- 
neous modulation in the magnitude of the exchange en- 
ergy which breaks lattice symmetries (as in spin-Peierls 
or 'plaquette' states), or 'fractionalization' leading to ex- 
citations characteristic of the deconfined phase of a Z2 
gauge theory ||, ||, ||, 0, [|, g|. However, general princi- 
ples which place constraints on the possible co-existence 
of such orders, on phases which can be separated by a 
continuous quantum phase transition, and on the topol- 
ogy of the phase diagram have so far not been delin- 
eated. Some issues along these lines have been addressed 
in recent studies |h], |ll), where the problems have typi- 
cally been related to strongly- coupled field theories whose 
properties are not reliably understood. 

It is clear that large scale numerical studies will be re- 
quired to address these subtle questions, but these have 
not been possible so far, because Monte Carlo simula- 
tions are impeded by the complex weights of the quan- 
tum spins. Exact diagonalization and series expansion 
studies yield only limited information which cannot eas- 
ily describe the vicinity of phase boundaries. One strat- 
egy to overcome these obstacles is to abandon study of 
the quantum spin model itself, and to focus instead on 
effective lattice models which can more easily explore the 
different limiting regimes of the phase diagram|l2|]. We 
will follow such an approach here. We will consider only 
quantum antiferromagnets with an easy-plane plane U(l) 



symmetry, as that makes them amenable to duality map- 
pings: the dual representation of our lattice model has 
only positive weights, allowing Monte Carlo exploration 
of the phase diagram on reasonably large lattices. This 
will allow us to present the first results on the phase di- 
agram of a model containing states whose orders extend 
over all those noted above. 

We will study an effective lattice model defined on 
a 3-dimensional cubic lattice of sites, j, representing a 
discretized spacetime with the three directions x,y,T, 
the last being imaginary time. Each site has an angu- 
lar degree of freedom, 6j, with (— e t9j ~ S+j, the 
5 = 1/2 spin raising operator on site j; the prefactor of 
e l9j is the sublattice staggering associated with the easy- 
plane Neel ordering. In addition, there is a Z2 gauge field, 
Sj,j+p, = ±1 residing on the links of the cubic lattice, 
where \x = x,y,r; this is an auxilliary field which arises 
in a number of distinct derivations of the effective model 
from the underlying quantum spin model ^ ^, p"3| , 
and including it allows for an especially simply descrip- 
tion of phases with fractionalization. The partition func- 
tion controlling these degrees of freedom is 

z= yi j n ° x p [ k e n s «+* 



{ S j.i+A =±1 } 



□ □ 



where A M is a discrete lattice derivative (A M /j = fj+p, — 
fj for any fj), and □ indicates all the elementary square 
plaquettes of the cubic lattice. A fundamental property 
of Z is its invariance under the Z2 gauge transforma- 
tion Sjj+ft — > WjSjj+firij+n, 6j — > 6j + 7r(l - rjj), where 
rjj = ±1 is an arbitrary field (we assume periodic bound- 
ary conditions along the r direction) . The model is char- 
acterized by two dimensionless coupling constants, g, and 
K, and we will explore its phase diagram in the g, K 
plane. The coupling g controls the propensity to Neel 
order in the easy-plane, with the small g region having 
magnetic order with (e tej ) ^ 0. The value of K is con- 
trolled most easily by ring-exchange terms in the under- 
lying quantum spin systemjn 



14, 15 . At K = 0, we 
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can independently sum over Sj.j+/i on each link, leading 
to weights dependent only on cos(A M 0j), which describes 
bosonic, S z = 1 quanta of e 3 hopping in spacetime; con- 
versely at large K (strong ring exchange), the Z 2 flux of 
the gauge field Sjj+p, is suppressed, and the 1/g term 
in describes the propagation of S z = 1/2 quanta, 
e i0j/2 ('gpinong'^ allowing the possibility of fractional- 
ized phases. The last term in (Q) is the crucial remnant 
of the Berry phases of the S = 1/2 spins of the quan- 
tum antifcrromagnet: it imposes quantization of half- 
integral spin on each lattice site. Notice that this term 
can make the weights in ([!]) negative: this is the central 
reason for the novelty and difficulty of the theory of two- 
dimensional quantum antiferromagnets. The theory Z 
was proposed as an effective model for S = 1/2 quantum 
antiferromagnets in Ref. [l^; it was also obtained Jl^| by 
taking the easy-plane limit of the effective lattice mod- 
els of quantum antiferromagnets with SU(2) symmetry 
proposed some time ago fl^| . Upon identifying the spin 
raising operator with a hard-core boson S+j ~ b\ : we 
can also consider Z as a model of a boson system with 
off-site interactions and plaquette exchange terms, at a 
mean density per site which is half- integral [^6| . 

Sedgewick et aZ.jlTj have recently studied the proper- 
ties of Z but without the Berry phase term: their physi- 
cal motivation was slightly different (the e l9j represented 
Cooper pair quanta in a superconductor) , for which this 
neglect could be justified under suitable conditions. For 
our purposes of using Z as a model of quantum antifer- 
romagnets, it is essential to include the Berry phase, and 
a description of its consequences is one of the purposes of 
this paper. However, because of the non-positive-definitc 
weights, we cannot perform a direct Monte Carlo sam- 
pling over 9j and Sjj+fr configurations, which was the 
method used by Sedgewick et q/.|l7|. Instead, we will 
now present a dual representation of Z which does have 
positive-definite weights. 

The duality mapping jllj proceeds by first replacing 
each cosine term in Z by a periodic Gaussian: 



the degrees of freedom are integers tj^ on the sites, j, of 



the dual cubic lattice; tanh Kd 



-2K 



is the coupling 



W(6) 



(2) 



where the m JjU are integers on the links of the direct 
lattice, and the prime indicates the constraint Tnj tli is 
even (odd) for Sj^+p, = 1). In this form Z can be 

subjected to a series of mappings which are described 
in Ref. [!(]: these are standard in the theory of duality 
of models with Z 2 and U(l) symmetry and will not be 
reproduced here. The end result is the partition function 

V J □ J 

(3) 

in which 6j and Sj^+p have been integrated out, and 



dual to K, the second sum on □ is over plaquettes on 
the dual lattice, and e^ u \ is the antisymmetric tensor. 
The fields ejj+p,, 0j,j+n are Ising variables on the links 
of the dual lattice: £jj+/i takes a set of fixed values such 
that Yia £ JJ+fi- i s "1 on ever y spatial plaquette and +1 on 
all other plaquettes (these values are linked to the Berry 
phase in (|l|)), while crj,j+p, measures the parity of £j M : 



1 - 2(£ JjU mod 2). 



(4) 



The physical meaning of the fields in (|3|) becomes clear 
upon relating them to observable properties of the anti- 
ferromagnet. The quantity (\ /2)e liv \^. v tj\ is associated 
with links of the direct lattice, and is the conserved 3- 
current of the spin-flip bosons foj introduced earlier: no- 
tice that this current can take half- integral values, rep- 
resenting the transfer of S z = 1/2 spinons in a singlet 
background. Also, the fj, = t component measures the 
boson density, or S z j. The first term in (|J) lies on the 
links of the dual lattice, each of which can be associ- 
ated with a plaquette of the direct lattice; for ji = x,y, 
each such plaquette is associated with a single bond of 
the antiferromagnet. In this manner we conclude that 
Sjj+fiO'jj+ii for [i = x, y is a measure of the exchange 
energy of the antiferromagnet on corresponding bonds in 
the y, x directions. 

We begin our study of Zd by considering various lim- 
iting regimes where it reduces to models studied earlier. 
(i) Kd = oo (K = 0): In this case crj,j+/i = and (||) 

becomes exactly equivalent to a dual current loop model 
of the Hubbard model of the fot bosons at half-integral 
filling considered by Otterlo et al. [Q. This also in 
keeping with the our earlier discussion that this limit of 
([!]) describes the facile motion of e l8j quanta. From the 
earlier work [|l8) we know that the boson ground state 
is superfluid (i.e. the spin system has easy-plane Neel 
order) for all g in such a model with only on-site interac- 
tions, (ii) Kd = (K = oo): The term involving crjj+p, 
disappears, and (||) becomes a model of half-boson cur- 
rent loops with on-site interactions. This is dual jl8| to 
the quantum XY model of half-angles, Oj/2, obtained in 
(Q) in this limit. This model has a superfluid phase at 
small g, and a fractionalized phase with freely propagat- 
ing and gapped S z = 1/2 spinons at large g. (Hi) g = 0: 
Both (Q) and (^) become trivial, and the ground state 
has superfluid (easy-plane Neel order), (iv) g = oo: The 
action (|) is that of a pure Z 2 gauge theory with a Berry 
phase term. In (^), this limit implies that <Jjj+ji = ijij+fi 
with $j = ±1 an Ising spin on the sites of the dual lattice; 
then (pf) becomes an 2+1 dimensional Ising model with 
every spatial plaquette frustrated, introduced in Ref. || as 
a model for the non-magnetic states of frustrated quan- 
tum antiferromagnets. The c,j quanta are the 'visons' of 
Ref. ^. This model has a fractionalized phase at small Kd 
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FIG. 1: Monte Carlo results for the phase diagram of Zd 
as function of g and Kd- Thin lines represent second order 
transitions, while the thick line appears to be first order. We 
present evidence below that the hatched region has co-existing 
Neel and bond order. As Kd — > oo, the region with superfluid 
order broadens and eventually extends to all g at Kd = oo. 




FIG. 2: Monte Carlo results for the Binder cumulant, rs, of 
the bond order parameter ^b at g = 14. The results indicate 
a second order transition from a phase with {^b) = at 
small Kd (where the value of tb is that associated with a 
Gaussian random field ^ b with zero mean), to a phase with 
(^b) 7^ at larger Kd (where tb = 1, the value expected 
by ^ b condenses). The crossing points indicate that this 
transition occurs at Kd ~ 0.35. 



(visons gapped) , and a confining, bond-ordered phase at 
large Kd (visons condensed). 

Combining all the limiting cases above, and with infor- 
mation gained from Monte Carlo simulations on (||) to be 
described below, we obtain the phase diagram shown in 
Fig [l| The simulations were performed by the Metropo- 
lis algorithm on cubic systems of size L x L x L, with 
periodic boundary conditions. 

Let us consider first the behavior as a function of 
Kd for large g. At g — oo, as mentioned above, 
there is a second order transition^) from a fraction- 
alized phase at small Kd to bond-ordered phase at 
large Kd- Our simulations show that this transition 
persists into the region with g large but finite. As 
before^], we define a complex bond-order parameter 

*s = Ej[(- 1 F" £ jj+^M+*+ i (- 1 ) j! ' e J,j+fi a 'j,j+y]; this 

order parameter is non-zero for both columnar and pla- 
quette orderings of the bonds, which are distinguished fief] 
by distinct values of arg(( , 5 b})- We merely tested 
for non-zero values (^b) and did not perform the 
more subtle tests required to determine ang^^S b))', this 
we did by computing the Binder cummulant tb = 
(I*s| 4 )/(I*b| 2 ) 2 - The results are shown in Fig |, and 
lead to the roughly horizontal phase boundary between 
the bond-ordered and fractionalized phases in Fig ^. 

Next we looked at the boundary between the Neel 
and bond-ordered phases as a function of g for differ- 
ent values of Kd- The simplest mean- field theory of Zd 
predicted Q that this should be a direct second-order 
transitions. Typical sets of our numerical results are 
shown in Figs [5] and ^. Apart from rs, we measured 
two additional physical quantities. The quantity Cy is 
the root-mean-square fluctuation in the action per unit 
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FIG. 3: Monte Carlo results for Kd = 1.0 as a function of g. 
The dashed lines indicate positions of proposed phase transi- 
tions. Notice that there is a small peak in Cv for the largest 
system sizes at the transition at the smaller value of g. Data 
for p s are not available for the largest system sizes because of 
long equilibration times in the winding number sector. 



4 



K H =0.5 




proportional to the system volume: this latter observa- 
tion suggests a first-order transition. 

This paper has described the phase diagram of a 
2+1 dimensional XY model coupled to a Z 2 gauge field 
with a Berry phase term. We reviewed arguments 
{§ [10[ |]1> H H asserting that this is an effective model 
for frustrated, easy-plane, S = 1/2, two-dimensional 
quantum antiferromagnets. We believe this is the first 
study of quantum phase transitions in the presence of 
Berry phase terms (which completely accounts for fluc- 
tuations of gauge fields), and which allows for the various 
orders that may be present in the ground state: magnetic, 
bond, and fractionalization. We found that fluctuations 
induced an intermediate regime with co-existence of mag- 
netic and bond orders in our U(l) symmetric model; a 
similar co-existence has also been argued recently for 
frustrated S = 1/2 models with SU(2) symmetry |u| |9). 

We thank T. Senthil and M. P. A. Fisher for very useful 
discussions. This research was supported by US NSF 
Grant DMR 0098226. 



FIG. 4: As in Fig | but for K d = 0.5. The inset shows the 
system volume dependence of maximum value of CV at the 
peak as a function of g. The proportionality is suggestive of 
a first order transition. Notice that the L dependence of the 
peak in CV is much stronger than that in Fig ^[ 



volume: if we interpret Zd in (g) as a classical statistical 
mechanical model in 3 dimensions, then Cy would be its 
specific heat (note, however, that CV is not the specific 



heat of the quantum antiferromagnet). It is an unbiased 



measure of the location of phase transitions and of their 
order. In addition we measured the stiffness, p s of the 
Neel (or superfluid) order: this was done, as in Ref. [t|, 
by generalizing the current loop model to its non-zero 
'winding number' sector, and measuring the fluctuations 
in winding number. 

T he weight of the evidence in Fig for Kd = 1.0 favors 
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two lsecond order transitions as a function of q, which are 



denoted by the dashed lines. At small g we are clearly 
in the superfluid phase with a non-zero p Sl and no bond 
order. Conversely at large g, bond order is present, but 
there is no superfluidity as p s is vanishingly small. How- 
ever the onset of superfluid order (as measured by p s ) 
clearly occurs at a point distinct from the onset of bond 
order (as measured by r^): there is an intermediate range 
where superfluid and bond order are both present. This 
is also supported by the measurements of Cy , which show 
pea ks at both transitions. 
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